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Abstract 

We show that non-frustrated and frustrated ladders in a magnetic field can be sys- 
tematically mapped onto an XXZ Heisenberg model in a longitudinal magnetic field 
in the limit where the rung coupling is the dominant one. This mapping is valid in 
the critical region where the magnetization goes from zero to saturation. It allows 
one to relate the properties of the critical phase {H}., H^, the critical exponents) to 
the exchange integrals and provide quantitative estimates of the frustration needed 
to create a plateau at half the saturation value for different models of frustration. 

PACS Nos : 75.10.jm 75.40.Cx 75.50.Ee 
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I. INTRODUCTION 



Intermediate between ID and 2D, ladders have been the subject of an impressive amount 
of work over the past few years [|I|. Thanks to an intensive experimental and theoret- 
ical effort, quite a lot is understood concerning the properties of S=l/2 ladders in a 
magnetic field. In particular, the magnetization starts to increase above a magnetic field 
and saturates above a magnetic field H^, and the phase realized for intermediate magnetic 
fields is believed to be a Luttinger liquid with gapless excitations and a power law decay of 
the correlation functions. As usual, it is difficult however starting from a microscopic de- 
scription in terms of exchange integrals to calculate the parameters of the low energy theory 
in the Luttinger liquid phase, and this description is to a certain extent phenomenological. 

In this paper, we show explicitly how this problem can be mapped onto the XXZ model 
in a longitudinal magnetic field if the rung coupling is the dominant one. This is actually the 
case of Cu2(C5lIi2N2)2Cl4, the ladder system on which most results under strong magnetic 
fields have been obtained so far Such a mapping can actually be performed for any 

type of coupling between the rungs, and we will study the following Hamiltonian (see Fig. 

1) 

N ^ ^ N 2 ^ ^ N ^ ^ 

1=1 1=1 a=l 1=1 

N ^ ^ N 2 

+ J2^Si^2Si+i,i- Sl^ (1) 

i=l i=l a=l 

In this expression, a (resp. i) is a chain (resp. rung) index, is the total number of 
rungs, and periodic boundary conditions along the chain direction are implicit. If all the 
couplings except are set to zero, the system is a collection of independent rungs. The 
states of a given rung are denoted by 15 >= (| ti> — I []>)/^/^^ \Ti >= \ TT>) l^o >= 
(I Ti> +1 iT>)/v^ and |T_i >= | |i>. In a configuration |cricr2 >, ai (resp. 0-2) refers to 
chain 1 (resp. 2). Their energies are E{S) = -3J±/A, E{Ti) = JJA - H, E{To) = J^/A 
and E{T^i) = J±/A + H. So upon increasing the magnetic field the groundstate of a given 
rung undergoes a transition between the singlet \S > and the triplet |Ti > at = J±, and 
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the total magnetization of the system jumps discontinuously from zero to saturation. 

If the other couphngs are non-zero but small, this abrupt transition is expected to broaden 
between and H^, — Hi being of the order of the largest of the couplings Ji, J2 and 
J2 • In this limit, the properties of the system for H] < H < H^ are best understood by 
splitting the Hamiltonian into two parts: 

N ^ ^ N 2 

'Ho = Jj_Yl Si,lSi,2 ^ -^c X! X! ^la^ 
1=1 1=1 a=l 

N 2 ^ ^ N ^ ^ N ^ ^ N 2 

H-l = X/ X/ ^i,aSi+l^a + ^2 X/ ^iA^i+i,2 + <^2 X/ ^i,2Si+l,l — {H — H^) ^ ^ Sf^^ (2) 
i=l a=l i=l i=l i=l a=l 

The groundstate of Tio is 2^ times degenerate since each rung can be in the state IS* > 
or |Ti >, and the first excited state has an energy equal to Jj_. Hi will lift the degeneracy 
in the groundstate manifold, leading to an effective Hamiltonian that can be derived by 
standard perturbation theory. Let us start by introducing pseudo-spin S=l/2 operators 
that act on the states IS* >« and |Ti of rung i according to 

^^^1-5 >.= -^\S >, <|Ti >,= ^|Ti >, 
a+\S >i= \Ti >i a+\Ti>i=0 

ar|5>,= cr^\Ti>,= \S>, (3) 

Then, to first order, and up to a constant, the effective Hamiltonian reads: 

N N 
i=l i=l 

The parameters of i^es are given by 

T T" 

- ^1 2 2 

"2 + 4 + 4 

H'^^ = H-H,-^--^--^ (5) 
^ 2 4 4 ^ ' 

The Hamiltonian of Eq.(^) is nothing but the XXZ model in a longitudinal magnetic 

field. This problem has been studied by several authors over the years, and most of the 



relevant information concerning the properties of the model is available in the literature 
0. In particular, the exponents of the spin-spin correlation functions have been obtained 
analytically when the model is integrable and numerically otherwise. To translate these 
results into the language of the original Hamiltonian of Eq.(|l|), one just has to express the 
original operators S^^, and S^^ in terms of the pseudo-spin operators. This can be done 
by inspection, and the results are: 



1 , 1 



(J, 



s/2 ' V2 

Sh = \i<yt + \) St^ = \{al + \) (6) 

A detailed discussion of the correlation functions measured in NMR experiments can be 
found in Ref. 0. We now discuss the implications of this mapping in different cases. 



II. THE REGULAR LADDER: = 4' = 

In that case J'^ = Ji is twice as large as J|^: We are in the XY universality class of the 
anisotropic Heisenberg model. The system is gapless, and it behaves as a Luttinger liquid. 
Having explicit expressions of the coupling constants in terms of the microscopic parameters, 
one can calculate everything in terms of these parameters. For instance, we can express 
and in terms of Jj_ and Ji. This is most easily done by first performing a Jordan- Wigner 
transformation to map the problem onto a problem of interacting, spinless fermions: 

N N N 

T^SF = i + h.c.) + VJ2 niUi+i - fi^rn (7) 

i i i 

The parameters of this Hamiltonian are given in terms of those of Eq.(^) by t = J!^y/2, 
V = Jf^ and yU = H'^^ + J^^. Hi corresponds to the chemical potential at which the band 
of spinless fermions starts to fill up. In that limit the repulsion term is irrelevant because 
the density of spinless fermions vanishes, so that the chemical potential corresponding to 
Hi is given by /i = — 2t. This leads to the result Hi = J± — Ji- To estimate H^, one cannot 
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neglect the repulsion term because the band is completely filled. The simplest way to take 
it into account is to perform a particle-hole transformation on the Hamiltonian of Eq.(0): 
c| — > (ij. Up to a constant, the new Hamiltonian reads 

N N N 

Hhole = -t E(4rf.+i + h.c.) + VY: nfnl, - fi^E^f (8) 

i i i 

where the hole chemical potential /i/i is given by = — + 2V. In terms of holes, corre- 
sponds to the chemical potential where the band starts to fill up, and one can again neglect 
the repulsion term. Note however that this is not equivalent to neglecting the repulsion in 
Eq.(^ since V appears in the expression of fih- The chemical potential corresponding to 

is thus given by fih = — 2t, leading to = J± + 2Ji. These expressions of and 
agree with those of Ref. obtained along different lines, and they compare well with the 
experimental values for Cu2(C5Hi2N2)2Cl4 0. 

The same argument actually apply if J2 and are not equal to zero. To first order, 
is unaffected, and is given by = J_l — Ji + (Jg + J2 



III. THE FRUSTRATED LADDER: = = J2 

When J2 7^ 0, the effective Hamiltonian is in the universality class of the XY model only 
if J2 is not too large: There is a transition to the Ising universality class when JJ^ = Jf^, 
i.e. J2 = ^1/3 to first order. In terms of spinless fermions, the Ising limit means that V 



is large enough to make the half-filled system insulating JTOl. The chemical potential as a 
function of the band-filling will then have a jump. In the original spin language, this implies 
that there will be a plateau in the magnetization at half the saturated value as a function 
of magnetic field. Note that similar conclusions have been obtained along different lines 



by several authors |Tl|- p!^] concerning the case J2 = (see next section). In the plateau 
region, there is an order parameter corresponding to alternating singlets and triplets on 
neighbouring rungs. For J2 = Ji, the effective Hamiltonian becomes purely Ising: Jf^ = Ji, 
J^y = 0. This result, clearly valid up to first order after Eq.(^, is actually exact including 
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all order corrections. The simplest way to understand this is to realize that a singlet on a 
given rung is completely decoupled from the rest because all the exchange integrals starting 
from this singlet belong to a pair of equal exchange integrals connecting a spin to both ends 
of the singlet. So the perturbation cannot couple to the singlets, and the XY exchange 
integral must vanish. Besides, this argument shows that the state with alternating singlets 
and triplets is an eigenstate of the Hamiltonian. It is easy to prove that it is the groundstate 
for H between and H^, which are given by = J± and = J± + 2Ji in the present 
case. So the plateau will extend over all the intermediate region between zero and saturated 
magnetization. 



IV. THE ZIGZAG LADDER: = 

This case is very similar to the previous one. There will be a transition to an Ising 
phase when is large enough, a conclusion already reached by other authors using different 



arguments |Tl|,|T2|JT^ . To first order the critical value is given by = 2Ji/3. The only 
difference is that there is no pure Ising phase in that case since a singlet is only decoupled 
from neighbouring singlets, and not from neighbouring triplets. The reason for mentioning 
this particular case of frustration is that it corresponds in principle to the physical situation 
realized in Cu2(C5Hi2N2)2Cl4. Our estimate of the critical value of to enter the Ising 
phase can be used as an upper bound to this exchange integral in Cu2(C5Hi2N2)2Cl4 since 
no plateau at half the saturated value has been reported. With Ji = 2.4 K, this means that 
J2 cannot exceed 1.6 K. Although there is some discussion in the literature as to what the 
actual value of this parameter is, all the estimates reported so far appear to be smaller than 
this upper bound. 



V. CONCLUSION 



In conclusion, we have shown that a strong coupling approach starting from the limit of 
strong rungs provides a simple and unifying picture of the very rich physics that appears 
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when ladders are put in a magnetic field. On one hand, it gives a simple explanation of how 
the Luttinger liquid physics emerges in the intermediate phase of unfrustrated ladders. On 
the other hand, this approach naturally leads to the presence of a plateau in the intermediate 
phase at half the saturation value when the coupling between the rungs is strongly frustrated. 
In systems where is effectively the largest coupling, this calculation allows one to relate 
measurable quantities like H^, and the critical exponents of the spin-spin correlation 
functions to the exchange integrals. Reported values for and in Cu2(C5Hi2N2)2Cl4 
are well reproduced by this approach. It will be interesting to analyze the critical exponents 
along the same lines when experimental data are available. Finally, this approach provides 
quantitative estimates of the frustration needed to create a plateau at half the magnetization 
value for systems where the rung coupling is the largest one and should help in the search 
for systems exhibiting this remarkable property. 

I acknowledge useful discussions with C. Berthier, M. Horvatic, M.-H. Julien, L. Levy 
and D. Poilblanc. I am especially indebted to Dr. A. Furusaki for pointing out a mistake 
in an early version of this paper. After completion of this manuscript, I have been informed 
by L. Levy that his group has independently obtained the same effective Hamiltonian in the 
case of unfrustrated ladders ||TB| . 
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FIGURES 




FIG. 1. Sketch of the ladder considered in this paper. 



8 



REFERENCES 

[1] For a review, see E. Dagotto and T. M. Rice, Science 271, 618 (1996). 
[2] P. R. Hammar, D. H. Reich, J. Appl. Phys. 79, 5392 (1996). 
[3] G. Chaboussant et al, Phys. Rev. B 55, 3046 (1997). 
[4] G. Chaboussant et al, Phys. Rev. Lett. 79, 925 (1997). 
[5] G. Chaboussant et al, unpublished. 

[6] C. Hayward, D. Poilblanc and L. Levy, Phys. Rev. B 54, R12649 (1996). 

[7] R. Chitra and T. Giamarchi, Phys. Rev. B 55, 5816 (1997). 

[8] Z. Weihong, R. P. Singh and J. Oitmaa, Phys. Rev. B 55, 8052 (1997). 

[9] F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

[10] A. A.Ovchinnikov, it Zh. Eksp. Teor. Fiz. 64, 342 (1973) [Sov. Phys. JETP 37, 176 
(1973)]. 

[11] K. Totsuka, to appear in Phys. Rev. B. 
[12] K. Totsuka, J. Magn. Magn. Mat. 177-181 (1998). 
[13] D. Cabra, A. Honecker, P. Pujol, unpublished. 
[14] T. Tonegawa et al, to appear in Phys. Rev. B. 
[15] G. Chaboussant et al, unpublished. 



9 



